Comparing lattice Dirac operators with Random Matrix Theory* 

F. Farchioni, I. Hipt and C. B. Lang 

Institut fur Theoretische Physik, Universitat Graz, A-8010 Graz, Austria 

We study the eigenvalue spectrum of different lattice Dirac operators (staggered, fixed point, overlap) and 
discuss their dependence on the topological sectors. Although the model is 2D (the Schwinger model with 
massless fermions) our observations indicate possible problems in 4D applications. In particular misidcntification 
of the smallest eigenvalues due to non-identification of the topological sector may hinder successful comparison 
with Random Matrix Theory (RMT). 



1. INTRODUCTION 

In recent work we have been studying various 
aspects of the lattice Schwinger model . This 
model is a 2D U(l) gauge theory of photons and 
one or more fermion species. Of particular in- 
terest is the situation of massless fermions. In 
the quantized theory chiral symmetry is broken 
by the anomaly. The one flavor-model should ex- 
hibit a bosonic massive mode. 

For the non-perturbative lattice formulation 
chirality is a central issue. The Wilson Dirac 
operator explicitly breaks chiral symmetry. The 
Ginsparg- Wilson condition [Bj defines a class of 
lattice actions with minimal violation of chiral- 
ity. An explicit realization is Neuberger's overlap 
Dirac operator [Q. In another approach one at- 
tempts to construct so-called quantum perfect ac- 
tions, or fixed point actions (classically perfect ac- 
tions) [5j, also obeying the Ginsparg- Wilson con- 
dition [£]. 

In the Schwinger model framework we have 
been studying several of these suggestions. In 
J?| the (approximate) fixed point Dirac operator 
was explicitly constructed. It has a large num- 
ber of terms but has been shown to have excel- 
lent scaling properties for the boson bound state 
propagators. This is not the case for the Neu- 
berger operator ||; there scaling is not notice- 
ably improved over the Wilson operator. The 
overlap operator has eigenvalues distributed ex- 
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actly on a unit circle in the complex plane; for 
the (approximate) fixed point operator our study 
shows small (with smaller f3 = 1/g 2 increasing) 
deviations from exact circularity. In both cases 
we could identify chiral zero modes. Their oc- 
currence was strongly correlated to the geomet- 
ric topological charge of the gauge configuration 
^geo = y~ J2 X Imln[/i2(a;) (henceforth called v 
for brevity) with a rapidly improving agreement 
with the Atiyah-Singer Index Theorem (inter- 
preted on the lattice) towards the continuum 
limit. 

Studying the spectra of the Dirac operators 
suggests comparison with Random Matrix The- 
ory (RMT)||. There the spectrum is separated 
in a fluctuation part and a smooth background. 
Exact zero modes are disregarded. The fluctu- 
ation part, determined in terms of the so-called 
unfolded variable (with average spectral spacing 
normalized to 1), is conjectured to follow predic- 
tions lying in one of three universality classes. 
For chiral Dirac operators these are denoted by 
chUE, chOE and chSE (chiral unitary, orthogonal 
or symplectic ensemble, respectively) |§ . Various 
observables have been studied in this theoretical 
context. Comparison of actual data should verify 
the conjecture and allows one to separate the uni- 
versal features from non-universal ones. In par- 
ticular it should be possible to determine in this 
way the chiral condensate. 

On one hand the limiting value of the density 
for small eigenvalues and large volume, 
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provides such an estimate due to the Banks- 
Casher relation. This information is contained 
in the smooth average (background) of the spec- 
tral distribution. However, also the fluctuating 
part, in particular the distribution for the small- 
est eigenvalue P(X m - m ) contains this observable: 
Its scaling properties with V are given by unique 
functions of a scaling variable z = AFE, depend- 
ing on the corresponding universality class. Usu- 
ally this is the most reliable approach to deter- 
mine S, which then serves as an estimate for the 
infinite volume value of the condensate in the chi- 
ral limit. This method does not involve unfolding, 
averaging or extrapolation. 

Here we concentrate on our results for the stag- 
gered Dirac operator. It is anti-hermitian and (for 
m — 0) its spectrum is located on the imaginary 
axis, but it has no exact zero modes. RMT pre- 
dictions for the staggered action and the trivial 
topological sector have been confirmed also in 4D 
lattice studies jl0|. Here we emphasize, however, 
the role of non-zero topological charge. 

2. METHOD AND RESULTS 

In our study we construct sequences of (5000- 
10000) uncorrelated quenched gauge configura- 
tions for several lattices sizes (16 2 , 24 2 , 32 2 ) and 
values of (3 (2, 4, 6). For these sets we then deter- 
mine the various Dirac operators and study their 
spectral distribution. This way we can compare 
directly the effect of identical sets of gauge config- 
uration on the fermionic action. In |Q we discuss 
our results for the Neuberger- and the fixed point 
operator. Since these spectra lie on or close to 
a circle in the complex plane, one has to project 
them to the (tangential) imaginary axis. We find 
that they exhibit the universal properties of the 
(expected) chUE-class, unless the physical lattice 
volume is too small. 

In Fig. |l| we demonstrate the relevance of topo- 
logical modes. The e.v. distribution density is 
first shown without distinguishing between differ- 
ent v and we notice a pronounced peak at small 
eigenvalues. Splitting the contributions accord- 
ing to \v\ = and 1 we observe, that the peak is 
due to the non-trivial sectors v ^ 0. The trivial 
sector has a behavior typical for the shapes pre- 
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Figure 1. Distribution density of e.v.s of the stag- 
gered Dirac operator, (out of all 10000 configu- 
rations 29% have v = and 45% have \v\ = 1.) 

dieted from chRMT. For larger j3 and V the peak 
becomes more pronounced, justifying the hypoth- 
esis that it represents the "would-be" zero modes. 

Since RMT discusses the distribution excluding 
exact zero-modes we expect problems whenever 
one is in a situation without possibility to sepa- 
rate topological sectors (upper-most figure in Fig. 
|l|), if one then tries to represent the distribution 
for the smallest observed eigenvalue by chRMT 
functions. This is demonstrated in Fig. || where 
we plot the histograms for the smallest and the 
2nd smallest (shaded histogram) eigenvalues in 
the \v\ = 1 sector. For small 0, strong coupling, 
the histogram for the smallest e.v. behaves like 
the v = sector prediction. For large (3 the 2nd 
smallest e.v. follows a distribution expected for 
the smallest e.v. in the \v\ = 1 sector. 

The level spacing distribution (determined in 
the unfolded variable) clearly has chUE (Wigner 




Figure 2. Distribution for the smallest and the 2nd smallest (shaded histogram) eigenvalues in the v = 1 
sector, quenched (Nf — 0). The curves give the chUE predictions for the smallest e.v. in the v — sector 
(dashed line) in the v — 1 sector (full line). 
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Figure 3. Level spacing distribution (for the un- 
folded variable s) for L = 16, j3 = 4, compared 
to chRMT predictions for chUE (full), chOE 
(dashed) and chSE (dotted). 

surmise) shape (Fig. |J) for all sizes and (3. 

Having all eigenvalues we can of course calcu- 
late the fermion determinant for every gauge con- 
figuration and include dynamic fermions by ex- 
plicit multiplication. These "unquenched" results 
will be presented elsewhere. 
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